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ABSTRACT 

This paper provides testable predictions about starlight polarizations to con- 
strain the geometry of the Galactic magnetic field, in particular the nature of the 
poloidal component. Galactic dynamo simulations and Galactic dust distribu- 
tions from the literature are combined with a Stokes radiative transfer model to 
predict the observed polarizations and position angles of near-infrared starlight, 
assuming the light is polarized by aligned anisotropic dust grains. SO and AO 
magnetic field models and the role of magnetic pitch angle are all examined. All- 
sky predictions are made, and particular directions are identified as providing 
diagnostic power for discriminating among the models. Cumulative distribution 
functions of the normalized degree of polarization and plots of polarization posi- 
tion angle vs. Galactic latitude are proposed as tools for testing models against 
observations. 

Subject headings: galaxies: magnetic fields - ISM: magnetic fields - polarization - 
radiative transfer 
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Introduction 



A gap exists between theories of the large-scale Galactic magnetic field and observations 
that probe the same field. To bridge this gap, predictions from existing Galactic dynamo 
models that can be tested with current observational techniques, in particular polarization 
of background starlight, must be generated. 

Extensive theoretical work over the past fifty years est ablished the origin and expected 



nature of galaxy-scale dynamos. 



Kulsrud fc Zweibell (120081 ) provide a comprehensive review 



of cosmic magnetic field research and highlight many of the existing questions related to 
galactic dynamos. A leading th eory for magnetic field generation in the Galaxy is the 



mean-field "a-fi" disk dynamo (JMoffatt 



1978 : 



Ruzmaikin et al 



19881). T his fundament al 



idea has led to many numeric a l simulations of Galactic magnetic fields (lElstner et al. 



1992 



Brandenburg et al. 



Hanasz et al. 



2009 



1992 



Moss et al. 



1993 



Ferriere Sz Schmitt 



2000 



Kleeorin et al. 



2002 



2003 



20101 ). As these ideas have developed, they became more 



sophisticated and added much to the understanding of the G alactic magnetic field. Recent 



advances include understan ding the role of magnetic helicity (JBrandenburg &: Subramanian 



2005 



2008 



Shukurov et al 



Moss et al. 



Hanasz et al. 



20091). 



20061), the interplay between the disk and halo fields (jMoss fc Sokoloff 



20101) , and the first global MHD simulation of a cosmic ray-driven dynamo 



The presence of magnetic fields in the Galaxy has be en directly observed via several 



obser vational methods: Zeeman spli tting of spec t ral lines (Verschuur 



19931 ). F araday rotation of pu lsars (ISmith 



1968 



Manchester 



sources (ICooper &: Price 



Wielebinski fc Shakeshaft 



9621) . and synchrotron emission (IWesterhout et al 



1962 



Yusef-Zadeh et al. 



1984 



1968; 



Crutcher et al. 



1974r) and extragalact ic 



1962 



19891 ) . Magnetic fields in the 



Hiltnei 



Galapcy have also been indirectly inferred from polarization of background starli ght ( 

19491 ; lHalllll949l ; iMathewson Sz Fordlll970l ) and polarized thermal dust emission fiHildebrand 
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19881 ). Zeeman splitting and polarized dust emission only probe high-density regions 
in interstellar clouds and therefore only weakly contribute to our understanding of the 
large-scale structure and origin of the Galactic magnetic field. Polarized synchrotro n 



emission has been used to study resolved magnetic fields in other galaxies (JBeck et al. 



1996 



and references there in) and in the Milky Wa y, in particular the vertical magnetic fields in 



the Galactic center (jYusef-Zadeh et al. 



198J). 



Much of the work attempting to understand large-scale magnetic fields in our Galaxy 
has used Faraday rotation of polarized emission from pulsars and extragalactic sources. 



Attempts have been made to fit 
the observed rotation measures ( 



Nota fc Katgertl 



2010 



empirical magnetic fie 



Rand fc Kulkarni 



Van Eck et al. 



1989 



models for the Galactic disk to 



Men et al. 



2008 



Sun et al. 



2008 



20111 ). These attempts to measure basic properties 
of the Galactic magnetic field using rotation measures have proved to be challenging, 



particular 



Han et al 



y m 



ight of arguments for a n antisymmetric (oddl (lAndreasvan fc Makarov 



19971 ). or symmetric (even) f lFrick et al. 



2001 



Sun et al 



in the disk, and for fields that cannot be so simply classified (JMao et al. 



2008h 



1988 



axisy mmetric field 



20101). Additional 



evidence suggests and that the observed rotation measure dist ribution may be strongly 



affected by nearby, small-scale structures (IWoUeben et al. 



2010h . 



Polarization of background st arlight has p reviously been used to measure the large-scale 



orientation of the magnetic field (JHeiles 



(defined as eiTctaja\Br/ Bj) of the field (JHeiles 



1976 



, to e stimate the curvature and pitch angle 



galaxies (JHough 



19961 ). and to study magnetic fields in other 



19961 ). However, it has mostly been applied to the s tudy of magnetic 



fields in and around star forming regions and individual c l ouds (e.g.. 



Goodman et al. 



1995 : 



Arce et al. 



1998; 



Tamura et al. 



2006 



2007 



Sato et a. 



Kandori et al. 



1985 



20071). By 



comparison, near-infrared (NIR) polarimetry of background starlight is a relatively untapped 
resource for the study of the large-scale Galactic magnetic field. Optical polarimetry has 



-5- 



typically been restricted to probing magnetic fields within ~1 kpc fJFosalba et al.l |2002| ) . 
limiting its ability to probe the large-scale structure of the magnetic field. However, NIR 
polarimetry is less extincted by dust and can probe the magnetic field along longer lines 
of sight. The regular component of the Galactic magnetic field may be better measured 



by these longer s i 



Rand fc Kulkarni 



ght. 



1989 : 



i nes (longer than the s cale length of the random component; 



Ohno fc Shibata 



19931 ) as the random magnetic components will 



tend to null averages. By studying the ma gnetic field in the outer Galaxy, where t 



le dust 



density is lower than in the inner Galaxy ( ISpergel et al. 



1996 : 



Drimmel fc Sperge 



and where there are fewer supernovae to create magnetic field disturbances (JHeiles 



2001 



19871), 



NIR polarimetry may be used to probe the quiescent, regular Galactic magnetic field to 
multi-kpc distances. 

Starlight polarimetry is insensitive to the directional polarity of magnetic fields, 
being sensitive only to their orientations. Therefore, starlight polarimetry cannot directly 
address questions about the nature of magnetic reversals seen in rotation measure studies. 
However, starlight polarimetry may shed light on the nature of the poloidal c omponent of 



the Galactic magnetic field, complementing previous rotation measure studies (JTaylor et al. 



2009 



Mao et al. 



20101 ). By observing the field projected onto the plane of the sky, this 



method is sensitive to the Galaxy's poloidal magnetic field component and allows the nature 
of that field to be measured. 

To advance the understanding of the Galactic magnetic field, strong constraints on 
theory are needed. In this paper, physics-based models of the Galactic dynamo, combined 
with models of the Galactic dust distribution, are used to make predictions for starlight 
polarizations that can be tested with current and emerging observational techniques. This 
will constrain the models, and therefore the physics, producing the Galactic magnetic 
field. The NIR stellar polarization mechanism and a Stokes radiative transfer model are 
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outlined in §2. Section 3 presents the results of the NIR polarization simulations using the 
models described in §2. In §4, future comparison of these predictions with NIR polarimetry 
measurements and the constraints that can be placed on Galactic dynamo models are 
discussed. This work is summarized in §5. 

2. Simulating NIR Polarimetry 

This section describes the assumptions used in simulating NIR stellar polarizations. 
The polarization mechanism is briefly discussed, as the focus of this paper is more about the 
observables and less so about the physical mechanisms causing polarization. The specific 
magnetic field and dust models are listed, followed by a brief discussion of the properties of 
the simulated stellar populations. Finally, the Stokes radiative transfer model is presented. 

2.1. Polarization Mechanism 



Polarization of background starlight is thought to arise from aniso tropic dust grai ns 



aligned in magnetic fields (for a comprehensive history of this idea see 



Theories for alignment n aechanisms include paramagnetic dissipat ion ( 



1951 



Purcell 



1975 



1979|). mechanical (ga s streaming) alignment 



agnetic dissipation (Jones fc Spitzer 



1967), and radiative torques 



Lazarian 



20031). 



Davis fc Greenstein 



Gold 



1952|). superparam- 



Dolginov fc Mytrophanov 



1976 



Draine fc Weingartnei 



1996 



1997 



Lazarian fc Hoang 



20071 ). With the exception of 



mechanical alignment, the prediction is that spinning dust grains will preferentially align 
with their long axes perpendicular to the local magnetic field. Unpolarized background 
starlight passing through a medium with aligned grains will experience a larger extinction 
cross-section perpendicular to the magnetic field than parallel to it, so that the transmitted 
light becomes weakly linearly polarized parallel to the direction of the magnetic field as 
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projected onto the plane of the skjl^ 



2.2. Magnetic & Dust Models 



Two broad classes of Galactic magnetic fields need to be considered: axisymmetric 
(ASS) an d bisymmetric (BSS), thoug h higher-order symmetries and combinations are 



possible ( IWielebinski fc Krause 



19931 ) ■ ASS and BSS fields are d ifferentiated by their 



symmetry with respect to the ro tation axis 



been observed in other galaxies (JBeck et al. 



Krause et al 



1989bl ) and both classes have 



19961 ) ■ Galactic magnetic fields are further 



distinguished by thei r disk symmetry, eit her symmetric ( even) or antisyrn metric (odd) 



across the midplane ( iKrause et al. 



1989b; see Table 1 in 



Beck et al. 



19961 ). Here, only 



ASS magnetic field models are considered; these are classified as AO (odd) or SO (even). 



This work wil 



also consider d isk-even, halo-od d (DEHO) magn etic fields which have 



observational ( ISun et al.ll2008l ) and theoretical ( JMoss et al. 



20101) foundations. 



The polarizing mechanism described 



lere is sensi tive not to the directional polarity of 



the magnetic field, but to its orientation (JHeiles 



19961 ). Consider the two simple magnetic 



field geometries in Fig. [T] one in which the magnetic field inside a spiral arm has the same 
directional polarity as the magnetic field outside the arm (ASS magnetic field), and one in 
which the magnetic field inside the arm has the opposite directional polarity of the rest of 
the Galaxy (BSS magnetic field). A distant star whose light passes through this spiral arm 
to an observer would gain identical polarizations in both cases because the orientation of 
the fields is identical even though the directional polarities may differ. Ergo, background 
starlight polarimetry is insensitive to magnetic field reversals. 



^Throughout this work, polarization position angles will be measured in the Galactic 
coordinate system. 



The predicted NIR polarizations were generated using simple magnetic fields 
configurations: axisymmeteric analytic magnetic fields; and axisymmetric dynamo- 
generated configurations from the literatu r e, in particul ar the dynamo simulation s of 



Ferriere &: Schmitt 



2000) and 



MossetaL 



mm . From 



Ferriere fc Schmitt 



(120001 ). the 



numerical outputs of three SO and three AO magnetic field models were obtained (D. 
Schmitt, private communication). These magnet ic field models ha d magnetic pitch angles 



Moss et al 



teOlOl ). the numerical outputs 



from 7.0 — 10.6° and are listed in Table [TJ From 
of three DEHO magnetic field models (513b, 511b, and 527b) were obtained (D. Moss, 
private communication), corresponding to the three models shown in Fig. 4 of that paper. 
This current work focuses on using these dynamo simulations as physically-motivated 
representatives of AO, SO, and DEHO Galactic magnetic fields. Three analytic, axisymmetric 
models (these can also be classified as SO models) were also considered: one consisting of 
purely azimuthal magnetic fields (5,. = 3^ = 0) and two possessing azimuthal and radial 



fields (but no B^ component) wit 



1 magnetic pitch a ngles mat ching those found for recen t 



models in the literature of 11.5° ( iBrown et al. 



20071) and 24° ( IRuiz-Granados et al.lboioh . 



Since polarization of background starlight is not sensitive to the directional polarity of 
magnetic fields the simulations used do not include magnetic reversals. 



Two different dust models (ISpergel et al. 



1996 



Drimmel Sz Spergel 



200l|) of the 



Galactic dust distribution were used, both derived from NIR and FIR emission measured 
by the COBE/DIRBE instrument for Galactic latitudes |6| < 30°. Since the observations 
simulated here extend to \b\ > 30°, the exponential forms utilized in the dust models 
were assumed to ad equately describe the Galactic dust distribution far from the midplane. 



Spergel et al. 



values. 



19961) imposed axis ymmetry, with fixed disk scale length and scale height 



Drimmel fc Spergel! (120011 ) invoked a more sophisticated model which included an 



axisymmetric disk with a linear flair, the local arm, four spiral arms, and a Galactic warp. 
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2.3. Stellar Population Sample 

Observed background starlight polarizations contain contributions from all aligned 
dust along the line of sight to the stars, so stars at different distances may sample 
different polarizing media and exhibit different polarization properties. Therefore, to 
accurately simulate real observations, the background stellar sources used here must be 



reasonabl y distributee 



mode 



of 



Robin et al. 



( JReyle et al. 



along and across Galactic lines of sight. The stellar population 



(120031 ) has been shown to agree with NIR star counts from 2MASS 



20091 ). so this model was used to generate 10' x 10' simulated star fields across 



the entire sky. Monte-Carlo selections of stars were drawn from extinction-free stellar 
probability distributions. Stellar H-band (1.6/im) photometry for each star was calculated 
using Johnson-Cousins filter transmission curves and we re extincted as suming a standard 



diffuse Galactic extinction of 1.8 mag kpc ^ at V-band (jWhittet 



20031 ) in the disk. The 



stellar distribution model returned stellar types, atmosphere parameters, mass, distances, 
Johnson- Cousins apparent magnitudes, apparent colors, and visual extinctions. Using these 
stellar samples, starlight polarization predictions were made for the entire sky, spaced 
every five degrees in Galactic latitude and longitude. To simulate realistic observational 
constraints, stars dimmer than H = 14 from the Monte-Carlo stellar population samples 
were excluded. 

2.4. Radiative Transfer Model 

The propagation of light from an unpolarized background star to an observer through 
aligned dust grains imposes weak linear polarization. To calculate the polarization 
observables P (the degree of polarization) and PA (the Galactic position angle of 



^http://model. obs-besancon.fr/ 
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polarization, measured toward increasing Galactic longitude from the Galactic North pole), 
unpolarized light was modeled as leaving each simulated star and propagating to the 



observer. As the starlig 



were calculated (Martin 



it pro pagates, the changes in the Stokes I, U, and Q parameters 



1974 ). All simulations remained in the optically thin regime with 
no contribution from polarized emission from any dust, as is appropriate for NIR wavelength 
observations of the Galactic disk. This assumption may not apply within 1° of the Galactic 
plane where the dust optical depth can become large. 

To illustrate aspects of the radiative transfer model, consider two limiting cases: a 
line of sight through two distinct polarizing regions whose magnetic fields are aligned with 
the same orientation, and a line of sight through two polarizing regions whose magnetic 
field orientation on the sky are aligned perpendicular to each other. In the first case, as 
the initially unpolarized background starlight passes through the first region a weak linear 
polarization (along the orientation of the magnetic field) is imposed. The light passes 
through the second polarizing region which further polarizes the background starlight 
along the same orientation. In the second case, the unpolarized background starlight again 
gains polarization from the first polarizing region. However, the second region depolarizes 
the starlight and the observer sees much weaker or unpolarized starlight. Because of the 
possible depolarizing effect of multiple polarizing regions, a single star may be a poor probe 
of the full nature of the magnetic field along a line of sight. Therefore, this work uses 
ensembles of stars (distributed with distance) to characterize each simulated sky direction. 

The Sun is assumed to be located 8.5 kpc from the Galactic center and in the Galactic 
midplane. For each simulated star, the entire line-of-sight from the star to the Sun was 
divided into 10 parsec blocks. This size was chosen to be smaller than the typical variations 
found in any of the magnetic and dust models. At each block, the local dust density was 
calculated, the magnetic field direction from the models (as seen from the Sun's position) 
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was projected onto the sky, and the changes in Stokes parameters were cal culated by 



integrating through that block. The changes in the Stokes parameters, from 



Lee&Draine 



( Il985l ). are described by the equations: 



d{Q/i) 2 0/2 iQ±-Q\ 

— e p IT a cos Zip cos 7 



ds 

-^epTca sin2%p cos -^ — —^ (2) 



ds 



di ^ 2 JQw + Q^) .^^ 



where e is the polarizing efficiency, p is the local dust density, a is the effective radius of the 
dust grains (assumed to be 0.1/im 'astronomical silicate'), ip is the direction of the magnetic 
field projected onto the plane of the sky, 7 is the angle between the direction of propagation 
of the radiation and the local magnetic field direction, an d Q is the sca ttering efficiency 
of the long (perpendicular) and short (parallel) axes from iDraind ( 19851 ). The polarizing 



efficiency of the dust is assumed to be identical throughout the Galaxy. 

The degree of polarization and position angle (PA) can then be calculated: 

P = ^U^ + Q2 (4) 

PA= - arctan(f//Q) (5) 

The uncertainty, and possible variation, of the polarizing efficiency of dust in the Galaxy 
affect the calculation of the absolute degree of polarization. Instead, the distributions of 
modeled polarizations were post facto normalized by the maximum calculated polarization 
in order to facilitate comparison with observed NIR polarization distributions. No additional 
constraints have been placed on the simulated star parameters, and all stars with H < 14 
are used in the following analysis. 

Simulations for all 24 models were generated across the entire sky for five degree steps 
in Galactic latitude and longitude. Each of the 24 models consists of one magnetic field 
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model plus one Galactic dust model. The model designations and their associated magnetic 



field and dust models are listed in Tab 



fl2000h SO models: models 7-12 are the 



e m Model numbers 1-6 are the 



Ferriere fc Schmitt 



Ferriere fc Schmitt 



are the 



Moss et al. 



teoOOh AO models; models 13-18 



(120101 ) DEHO models; models 19-24 are the axisymmetric analytic 
models. For each model, the P and PA for each star in each Monte-Carlo stellar sample was 
calculated. The mean PA and its dispersion were calculated for one 10' x 10' field of view 
at each sky grid direction for each model. Cumulative distribution functions (CDFs) of the 
polarizations from each model were tabulated. The polarization CDFs are shown in Fig. [2] 
and the quartile values of the cumulative P/Pmax distributions are presented in Table [T] 
For example, for model 1, 50% of the PjPmax values are smaller than 0.0038. 



Results 



The variation of simulated starlight polarization across the sky, centered on the 
Galactic anti-center, for a subset of the SO, AO, DEHO, and axisymmetric analytic magnetic 
field models are shown in Figs. [3] through [261 Predictions towards the Galactic center, 
and in particular in the Galactic plane, are likely less reliable than the Outer Galaxy and 
higher Galactic latitudes because of the effects of star formation and supernovae which may 
interact with the large-scale magnetic field. 



The SO (Fig. [3]), DEHO (Fig. [T5]), and axisymmetric analytic (Fig. [22]) magnetic field 
models all show similar structures, with magnetic field orientations in the disk parallel to 
the Galactic plane and magnetic nulls near Galactic longitudes 90° and 270°, which are 
along the toroidal magnetic component. This is significantly different from the predictions 
of the AO models shown in Fig. IH] where the field is predominantly perpendicular to the 
Galactic plane and the magnetic nulls occur at high Galactic latitudes. 



-13- 



Toward the i = 0° and 180° directions, many of the model predictions become 
degenerate. In particular, the three SO, three DEHO, and 6 = 0° analytic models 
make identical predictions toward these directions at all Galactic latitudes. The other 
axisymmetric analytic models vary with Galactic latitude because of projection effects alone. 
The i = 90° and 270° directions are also degenerate in all models except the 6 = 11.5°, 24° 
models (Models 21-24). Toward these directions, the field lies predominantly along the line 
of sight, so the toroidal component has a small projection onto the sky and the (mostly 
vertical) poloidal field projection dominates. This is also reflected in the weak polarizations 
predicted for these directions. In the models with large pitch angles {6 = 11.5°, 24°; Models 
21-24), the magnetic nulls shift in Galactic longitude (see Fig. |23]), due to the change in the 
magnetic field geometry. The amplitude of the longitude shift depends directly on the pitch 
angle. The AO models are almost entirely perpendicular to the Galactic plane for \b\ < 45°. 
For |6| > 45°, the projection of the radial magnetic field component begins to align with the 
Galactic plane. Because of all these effects, the directions near i = 0°, 90°, 180°, and 270° 
are not suitable for distinguishing among the various magnetic field model geometries. 

The simulated data were extracted from all the models at one example Galactic 
longitude {£ = 150°). As shown in Fig. [271 the SO models are seen to all have similar shapes 
in the (PA, b) plane, especially when compared t o the very different AO models shown in 

(120101 ) in Fig. [291 ha ye similar shapes t o the SO 



Fig. 1281 The DEHO models o 



models of 



Ferriere fc Schmitt 



Moss et al 



(|2000[ ). This is expected since the 



Moss et al. 



fl2010[ ) DEHO 



models have even magnetic fields in the disk and odd magnetic fields in the halo. Both dust 
density models fall exponentially with Galactic height, so the modeled polarizations are 
more sensitive to fields in and near to the disk and less sensitive to far halo (high latitude) 
fields. The higher PA dispersions seen toward Galactic mid-latitudes in the DEHO models 
may be caused by partial sampling of the halo-odd field. The 6' = 0° models (Models 19 
and 20 in Fig. [30l) are qualitatively similar to the SO and DEHO models, however the other 
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axisymmetric analytic models (in Fig. |30|) show a flattening of their PA vs. b for the higher 
magnetic pitch angles. By 9 = 24° (Models 23 and 24 in Fig [501) . the distribution of PA as 
a function of Galactic latitude is nearly flat, with mean PAs around 90°. In this case, the 
degenerate region that exists near i = 180° for 6 = 0° has shifted towards i = 150° (as seen 
in Fig. |25]). 



The degeneracies typically seen towards i = 90° and 270° for the even models with 
small magnetic pitch angles (Models 1-6, 19, and 20) was broken for larger pitch angles 
(Models 21-24), thus the longitude locati on of t h e pola rization nulls may be used as a tool 



for measuring the magnetic pitch angle. 



Heiled (119961 ) has previously used the location 



of these polarization null points to estimate the Galactic magnetic pitch angle. Towards 
i = 90°, projections of the magnetic pitch angle causes a distinctive change in the shape 
of the PA vs. Galactic latitude plot. This may allow observations in this direction to put 
constraints on the magnetic pitch angle of the Milky Way. It is important to note however, 
that CDFs for this region of the sky were dominated by low polarizations, and that many 
of the predicted stellar polarizations may fall below realistic detection thresholds. 

The choice of dust models has a non-negligible effect. As an example, consider models 



13 and 14 in Fig. [29 



field from 



Moss et al. 



Drimmel &: Spergel 



3oth models use the same model stars , the same model m agnetic 



f|2010[) . but use different dust models (jSpergel et al. 



1996 



versus 



200 ll ). The RMS variation of the mean PA values across the sky 



between the two models is 1.09°, which is small or comparable to the mean PA dispersion 
for each model (0.91° for Model 13, and 1.77° for Model 14). The most significant difference 
between the two models is found in the intrinsic PA dispersions, shown as a ratio in Fig. 
3T] for the one Galactic longitude i = 150°. The differences seen must be due to the 
dust model details. This effect was seen along all simulated Galactic longitudes. The 



Drimmel &: Spergell (120011 ) dust model (used for Model 14) includes the local arm, spiral 
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arms, and a Galactic warp, which were not included in the ISpergel et al.l ( 119961 ) model 
(used for Model 13). The choice of dust model also affects the tabulated CDF vales for 
the different models, as listed in Table [TJ For any of the four classes of magnetic fields 



used, there are diffe rences among the three 



Spergel et al 



the 



Spergel et al 



Drimmel fc Spergell (J200ll ) models and the three 



|l996) mode ls. For the AO, SO, and axisjTiimetric analytic magnetic fields, the 



Drimmel fc Spergel 



Spergel et al 



Drimmel &: Spergell (120011) models produce slightly higher normalized polarizations than 



9961) m odels. DEHO magnetic fields show the opposite effect, with the 



( 120011 ) models producing slightly lower normalized polarizations than 



(Il996[ ) based models. 



All of these simulations are optically thin, and so wavelength i ndependent. However 



19751). 



wavelength is known to affect the polarizing efficiency of dust grains ( iSerkowski et al. 
Under the assumption that there is a uniform wavelength dependence of polarizing efficiency 
along the lines of sight, scaling between wavelengths only requires use of a multiplicative 
factor between the predicted, normalized degree of polarization distribution and the actual 
polarization distribution. These effects will not change the observed PA for a given star 
since, under these assumptions, the polarizing efficiency does not affect the PA calculation. 

In Fig. |2l the CDFs of the P values in each model are shown, normalized by the 
star with the highest predicted polarization. Each model was normalized independently. 
The star with the highest predicted polarization value is typically 20-40 times the average 
modeled polarization. In 75% of the models, the star with the highest polarization was 
a simulated G2 supergiant at 13.31 kpc, and in 25% of the models it was a simulated F8 
supergiant at 13.51 kpc (all of the models used the s ame stellar p o pulat ion samples). These 



Robin et al. 



(120031 ) stellar population 



were some of the most distant stars returned by the 

model. The long distances between these stars and observer generated larger degrees of 

starlight polarization. These bright stars are exceedingly rare, but may be detected in 
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magnitude-limited polarization surveys covering large fractions of the sky. Since these high 
polarization stars make up such a small fraction of the total population, the CDFs shown 
in Fig. |5]only span the lowest ten percent of the entire CDF. 

A discontinuity is seen in the CDFs at approximately the 93rd percentile, where 
distant, high polarization stars occur. There are also differ e nces to be noted between the 



odd-numbered models, which used the 



the even-numbered models, which use the 



Spergel et al 



Drimmel &: Spergell (|200lh dust distribution, and 



(119961) dust dis tribution. For a 



given magnetic field geometry, the models using the iSpergel et al. 



(119961 ) dust distribution 



show predicted polarizatio ns typically about four times larger than the models using the 



Drimmel fc Spergell (120011 ) dust distribution. This effect is strongest for low polarization 
stars and weaker for the stars with the highest polarizations. Normalization by the 
maximum predicted polarization does not mask the effect. It is caused by differences in 
the normalization constants used for the two dust distributions. The dust column was 
calculated for every sample star using the two different dust distributions. The average 
dust column ratio was 4.2, which is approximately the maximum value in a corresponding 
polarization ratio distribution. The interplay between the exact dust distribution details 
and projected magnetic field geometry creates some depolarization, which causes the 
polarization ratio to decrease for some stars. 

In addition to an unknown polarization scaling factor, actual observations will suffer 
from censoring at the low P end of the CDF, as stars fall below polarization detection 
thresholds. Stars above detection thresholds can be used to calculate an appropriate scaling 
factor to match observations to the current simulated predictions. In practice, the censoring 
of P below an observational pfolarization limit would manifest as a sharp cutoff at the low 
polarization end of the CDF. The low and high polarization ends of the CDF may also 
suffer from undersampling. The 25th, 50th, and 75th percentile CDF polarization values 
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should therefore best serve as appropriate locations to calculate this scaling factor, and 
these values are tabulated for each of the model CDFs and presented in Table [TJ 

In a magnitude-limited polarimetric survey, the observed polarization CDF of all stars 
(including non-detections) above the observational apparent magnitude limit would be 
constructed. A predicted degree of polarization CDF in the same region of the sky, at the 
same limiting magnitude, should be simulated. The ratios of the degrees of polarization at 
a given percentile value in the observed and simulated CDFs will give the normalization 
factor. This normalization factor can then be used to calibrate the predicted degrees of 
starlight polarization to actual observations. 

4. Discussion 

The simulations presented here attempt to bridge the gap between theory and 
observation by making testable predictions for starlight polarization observations based 
on models of the large-scale structure of the Galactic magnetic field. These predictions, 
when combined with NIR polarization observations, can be used to test the models for the 
large-scale magnetic field and their physical underpinnings. Polarization of background 
starlight is insensitive to magnetic reversals and cannot address the presence or location 
of reversals. However, these polarizations can allow testing for the existence and nature 
of the poloidal component of the Galactic magnetic field, since all SO models yield similar 
polarization predictions across the sky and all AO models yield similar predictions that are 
distinctly different from the SO models. The predicted longitude shift of the polarization 
null points with magnetic pitch angle, demonstrated in Models 19-24 (a small shift is seen 
all of the models, however only Models 19-24 systematically vary the pitch angle), may 
facilitate an observationa l test for the Galactic magnetic pitch angle in the disk, similar to 
the work of iHeilesI (119961 ). This test is fundamentally different from previous attempts to 
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use Faraday rotation of pulsars anc 



( JRand fc Kulkarni 



1989 



Men et al 



extragalactic sources to model the magnetic pitch angle 



2008 : 



Nota fc Katgert 



2010 



Van Eck et al. 



201l|) 



The Galaxy was assumed to have an axisyrnmetric magnetic field based on theory 



( jRuzmaikin et al. 



1988; 



observational evidence (Rand fc Kulkarni 



Ruiz-Granados et al 



M81, 



Krause et al. 



Brandenburg et a 



1990; 



m 



Moss fc Brandenburg 



19891 : iRand fc Lvne 



1994 



1992) and 



Sun et al. 



2008; 



20101) . but there is sorn e evidence in other galaxies (partic ularly in 



1989al : ISokoloff et al. 



19921 ) that bisymmetric fields may exist. 



Beck et al. 



( 1l996l ) indicate that many of the galaxies showing evidence for bisymmetric magnetic fields 



also show evidence for galaxy interactions, so we might not expect t his magnetic fi eld 



structure for the Milky Way (based on a lack of major merger events; 



Gilmore et al. 



2OO2I) 



and the axisymmetric assumption may hold. Higher order az imuthal symmetr ies might also 



be possible, but their amplitudes should be relatively small (jBeck et al. 



19961) 



The measurements needed to distinguish among the various model simulations 
are well suited to NIR stellar polarimetry. The polarization m echanism used in the 



simulations works from the NIR through near-UV waveleng ths (jSerkowski et al. 



Codina-Landaberry &: Magalhaes 



1976 



Whittet et al. 



1975 : 



19921 1. However, NIR light is less 



attenuated by interstellar dust and can probe magnetic fields along multi-kp c scales, while 



th e shorter wavelengths o nly probe within about 1 kpc (jFosalba et al 



by 



Serkowski et al. 



2OO2I ). As described 



(Il975[ ). the polarization signal in the NIR is weaker than the visible by 



a factor of four or more. Ho wever, NIR polarimetric observations at th is level are possible 



with recent instrumentation (JKandori et al 



2006 



Clemens et al 



20071 ) and should soon 



provide data able to test models of the large-scale structure of the Galactic magnetic field. 

The simplest comparison may be through collecting the polarization behavior exhibited 
through a set of samples of all Galactic latitudes at a single Galactic longitude, as shown 
in Figs. [27] through |30l To best understand the poloidal component of the Galactic 
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magnetic field, these observations should be made in the outer Galaxy longitude ranges 
oi i = 110 — 160° OT i = 200 — 250°. These ranges avoid the degenerate regions that 
reside near i = 0, 90, 180, and 270° with a 20° buffer zone in either direction to allow for 
possible wander caused by any reasonable magnetic pitch angles (if present). To measure 
this magnetic pitch angle, a sample of Galactic longitud es shou l d be observed at a single 



Galactic latitude, say b = 15°, similar to the approach of 



HeilesI (119961 ). 



5. Summary 

Predictions for the polarization of background starlight, based on physics-driven 
Galactic dynamo simulations and empirical Galactic dust distributions, are presented. 
A Stokes radiative transfer model was used to predict the observed stellar polarization 
properties from Monte Carlo generated realistic stellar distributions in the Galaxy. A range 
of sky directions are suggested as particularly diagnostic regions for differentiating among 
Galactic magnetic field model geometries. Samples across many Galactic longitudes at one 
Galactic latitude, however, may be more appropriate for characterizing the magnetic pitch 
angle of magnetic field patterns in the Galaxy. 

The CDF of observed starlight polarizations is proposed as a tool for calibrating these 
models. Table [1] presents quartile values, for each model, of the normalized degree of 
starlight polarization. These measures are easily obtained from observations of starlight 
polarizations and can be used to determine a scaling coefficient between the predicted 
degree of polarization and observed degree of polarization. 

The role of the Galactic magnetic field in interstellar dynamics and star formation 
is only beginning to be understood, but new technical advances will permit probing this 
mysterious component of our Galaxy. The study of the Galactic magnetic field has been 
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dominated by Faraday rotation measurements, but new tools, such as NIR polarization of 
background starlight, complement these studies. 

The author thanks D. Schmitt and D. Moss for providing electronic versions of their 
simulations, and K. Ferriere, T. J. Jones, D. Clemens, and A. Pinnick for helpful comments 
and discussions. The author would also like to thank the anonymous referee whose 
comments significantly improved this work. This work was partially supported by NSF 
grants AST 06-07500 and 09-07790 to Boston University, D. Clemens PI. 



Table 1. Model Designations 



Model Number 



Magnetic Field 



Magnetic Dust Model ^ / Pmax Percentiles [xlO^ 



Pitch Angle Reference** 25th 



50th 



75th 



SO reference run'' 

SO reference run'' 

SO reference with a quenching'' 

SO reference with a quenching'' 

SO with vacuum BC'' 

SO with vacuum BC'' 



8.4° 
8.4° 
9.1° 
9.1° 
9.1° 
9.1° 



DS2001 
SMB96 
DS2001 
SMB96 
DS2001 
SMB96 



13.6 
40.2 
16.8 
48.1 
16.8 
47.0 



38.3 
106 
47.6 
126 

47.8 
124 



151 
367 
178 
417 
189 
427 



9 

10 
11 
12 



AO reference run'' 

AO reference run'' 

AO reference with a quenching'' 

AO reference with a quenching'' 

AO with vacuum BC'' 

AO with vacuum BC'' 



10.6° 
10.6° 
7.0° 
7.0° 
7.3° 
7.3° 



DS2001 
SMB96 
DS2001 
SMB96 
DS2001 
SMB96 



21.1 
53.5 
21.3 
54.0 
21.3 
54.0 



57.8 
157 
57.7 
157 
58.1 
158 



103 

248 
104 
249 
104 
251 



13 
14 
15 
16 
17 



C^ind = 0, i?„,fta(o = 300== 8.1° DS2001 17.0 48.3 188 

C^„«nd = 0, i?„,fta(o = 300^ 8.1° SMB96 47.6 125 429 

C^„^nd == 100, i?„jiaio = 300^= 13.0° DS2001 17.4 49.8 193 

Cn,^nd = 100, RaJ^alo = 300^= 13.0° SMB96 48.2 129 439 

C^„^nd = 200, i?„,,iaio = 300^= 21.0° DS2001 15.9 44.6 170 



Table 1 — Continued 



Model Number 



Magnetic Field 



Magnetic Dust Model P/Pmax Percentiles [xlO''] 



Pitch Angle Reference'^ 25th 



50th 



75th 



18 



Cwind — 200, Rahalo — SOC^ 



21.0° 



SMB96 



45.1 



119 



398 



19 


Analytic Model, 6* = 0° 


0° 


DS2001 


10.6 


30.0 


112 


20 


Analytic Model, 6* = 0° 


0° 


SMB96 


33.7 


88.1 


298 


21 


Analytic Model, 9 = 11.5° 


11.5° 


DS2001 


11.0 


31.4 


120 


22 


Analytic Model, d = 11.5° 


11.5° 


SMB96 


35.2 


91.7 


311 


23 


Analytic Model, d = 24° 


24° 


DS2001 


12.7 


35.0 


125 


24 


Analytic Model, B = 24° 


24° 


SMB96 


40.5 


105 


339 



^DS2001 



Drimmel fc Snergell pOOlh : SMB96 



Spergel et al. 



|l99(Th . 



^Ferriere fc Schmittl ( 2000 ) 



'- Moss et al 



mm 



to 

to 
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Fig. 1. — Cartoon illustrating the insensitivity of polarization of background starlight to 
magnetic polarity reversals, seen in a top-down view of a galaxy. A star (star symbol) 
emits light that passes through a spiral arm (deliniated by dashed lines) to an observer (eye 
symbol) for the cases of an axisymmetric ASS (left) and bisymmetric BSS (right) magnetic 
field. The arrows indicate the directional polarity of the magnetic field in each region. The 
observer would measure identical linear starlight polarization properties for the star in both 



cases. 



-24- 



1.0 
0.5 
0.0 

:;= 1.0 

JD 
D 

-Q 0.5 
o 

Q_ 



(L) 
> 



O 

3 



o 



0.0 



1.0 



0.5 



0.0 



1 Model 1 


r 

/ Model 7 


' Model 13 


Model 19 


/ Model 2 


/ Model 8 


/ Model 14 


/ Model 20 


/ Model 3 


r 

/ Model 9 


1 Model 15 


Model 21 


/ Model 4 


/ Model 10 


/ Model 16 


/ Model 22 


1 Model 5 


/ Model 11 


' Model 17 


( 

Model 23 


/ Model 6 


/ Model 12 


/ Model 18 


/ Model 24 



0.00 0.05 0.10 



0.00 0.05 0.10 

P/P 



max 



Fig. 2. — Cumulative probability distributions (CDFs) of stellar polarization for all 24 
models, normalized by the maximum predicted degree of polarization, shown for cumulative 
probabilities of to 0.1. These CDFs can be used to calibrate the predicted degree of 
starlight polarization with real observations. 
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Fig. 3. — Predicted mean Galactic P and PA for 10' x 10' fields on a 5 x 5° grid for the SO Model 1. At each grid 
location, the uniform length vectors represent the PA of the observed polarization and the background color represents 
the degree of polarization. All of the SO magnetic field Models produce a similar shape. Figs. 4-8, 10-14, 16-22, 24, 
and 26 are available in the online version of the Journal. 
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Fig. 4.— Same as Fig. El but for SO-type Model 2. 
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Fig. 5. — Same as Fig. |3l but for SO-type Model 3. 
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Fig. 6. — Same as Fig. |3l but for SO-type Model 4. 
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Fig. 7. — Same as Fig. |3l but for SO-type Model 5. 
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Fig. 8. — Same as Fig. |3l but for SO-type Model 6. 
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Fig. 9.— Same as Fig. |3l but for the AO-type Model 7. 
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Fig. 10.— Same as Fig. [3l but for the AO-type Model 8. 
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Fig. 11.— Same as Fig. [3l but for the AO-type Model 9. 
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Fig. 12.— Same as Fig. [3l but for the AO-type Model 10. 
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Fig. 13.— Same as Fig. [3l but for the AO-type Model 11. 
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Fig. 14.— Same as Fig. [3l but for the AO-type Model 12. 



P /2 



max 



en 

CD 
X) 

<D 

X) 



o 
o 

o 
O 




00 



~\\\ \ \ \ u ; / / / //^ 

-NNW \ \ \ 11 ; / / / / /yyy^ 

.N\\\ \ \ n M / / / / / / y y^ 

~-v\\\\\\\i I I / \ I I / / / / ^ ^ 



ppiir^-~iii~-::^?^iiiP 





N\\\ \ \ U I / / / / / ^y-^-^' 

N\\\ \ \ I I I ; / 1 / / /y^^- 

■vsvww \ \ I I J / / I / / / yy^^ 
xx\\\\\\ \ I I / / I I / / / / y y. 



-90 ^^^^^^^^^^^^^^^M 




~-^-v\\\\\\\\ I 1 I ) ( I I I / / / / y y-^- — 



90 180 270 

Galactic Longitude [deg] 



360 



Fig. 15.— Same as Fig. [3l but for the DEHO Model 13. 
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Fig. 16.— Same as Fig. [3l but for the DEHO Model 14. 
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Fig. 17.— Same as Fig. [3l but for the DEHO Model 15. 
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Fig. 18.— Same as Fig. [3l but for the DEHO Model 16. 
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Fig. 19.— Same as Fig. [3l but for the DEHO Model 17. 
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Fig. 20.— Same as Fig. [3l but for the DEHO Model 18. 
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Fig. 21. — Same as Fig. [3l but for analytic axisyinmetric Model 19. 
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Fig. 22. — Same as Fig. [3l but for analytic axisyinmetric Model 20. 
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Fig. 23. — Same as Fig. [3l but for analytic axisyinmetric Model 21. 
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Fig. 24. — Same as Fig. [3l but for analytic axisyinmetric Model 22. 
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Fig. 25. — Same as Fig. [3l but for analytic axisyinmetric Model 23. 
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Fig. 26. — Same as Fig. [3l but for analytic axisyinmetric Model 24. 
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Fig. 27. — Predicted mean Ga. 



models, from 



Ferriere fc Schmitt 



actic polarization PA in 10' x 10' fields for simulated SO 



(I2OOOI ). at i = 150°, as a function of Galactic latitude. The 



error bars represent the ±5cr PA dispersions for all stars in each Galactic latitude field. All 
of the SO magnetic field models produce similar shapes. 
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ig. 28. — Similar to Fig. |271 but for AO magnetic field models from 



Ferriere &: Schmitt 



2000l ). Some data points appear twice because of the 180° ambiguity in PA. 
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Fig. 29. — Same as Fig. [271 but for DEHO models from 



Moss et al. 



torn . 
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Fig. 30. — Same as Fig. |271 but for analytic axisymmetric models with different pitch angles. 



5F 



-53- 



4 F 



n 



< : 



2 



< 

CL 

b 




-100 



-50 50 

Galactic Latitude [b] 



100 



Fig. 31. — Ratio of Model 13 and Model 14 PA dispersions as a function of Galactic latitude 
for a = 150°. 
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